IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION,  VOL.  49,  NO.  10,  OCTOBER  2001 


1411 


High-Resolution  Phase  Sampled  Interferometry 
Using  Symmetrical  Number  Systems 

Phillip  E.  Pace,  Senior  Member,  IEEE,  David  Wickersham,  David  C.  Jenn,  Senior  Member,  IEEE,  and  Nathan  S.  York 


Abstract — This  paper  identifies  a  new  phase  sampling  in¬ 
terferometer  approach  that  can  be  easily  incorporated  into  the 
established  techniques  to  provide  a  high  resolution,  small-baseline 
array  with  fewer  number  of  phase  sampling  comparators.  The 
approach  is  based  on  preprocessing  the  received  signal  using 
symmetrical  number  systems  (SNS).  Antennas  based  on  both 
an  optimum  symmetrical  number  system  (OSNS)  and  a  robust 
symmetrical  number  system  (RSNS)  are  Investigated.  The  SNS 
preprocessing  is  used  to  decompose  the  spatial  filtering  operation 
into  a  number  of  parallel  suboperations  (moduli)  that  are  of 
smaller  computational  complexity.  A  much  higher  direction 
finding  (DF)  spatial  resolution  is  achieved  after  the  N  different 
moduli  are  used  and  the  results  of  these  low  precision  subop¬ 
erations  are  recombined.  By  incorporating  the  OSNS  or  RSNS 
preprocessing  concept,  the  field  of  view  of  a  specific  configuration 
of  interferometers  and  phase  sampling  comparator  arrangements 
can  be  analyzed  exactly.  The  OSNS  gives  the  maximum  dynamic 
range  or  number  of  spatial  resolution  bins  while  the  RSNS 
reduces  considerably  the  number  of  possible  encoding  errors. 
Experimental  results  for  both  a  5-blt  OSNS  and  a  6-bit  RSNS 
array  are  compared.  The  errors  in  the  encoding  of  the  direction  of 
arrival  are  quantified  for  both  architectures. 

Index  Terms — Directional  finding,  electronic  warfare,  symmet¬ 
rical  number  system. 

I.  Introduction 

Direction  finding  (DF)  systems  find  applications 
in  many  disciplines  including  geolocation,  navigation 
and  targeting  for  weapon  systems.  They  are  also  used  for 
tracking  by  law  enforcement,  military  personnel,  and  wildlife 
management.  Radio  DF  systems  commonly  use  phase  sampled 
linear  interferometer  arrays  to  measure  the  direction  of  arrival 
of  a  single  incident  planar  waveform  in  order  to  estimate 
the  geographical  bearing  of  the  emitter.  Each  interferometer 
consists  of  two  receive  elements  separated  by  a  distance  d 
(baseline).  A  waveform  impinging  on  an  ideal  interferometer 
has  the  same  amplitude  at  both  element  locations.  Thus,  all 
of  the  information  about  the  direction  of  arrival  is  contained 
in  the  phase  of  the  signal  measured  by  the  different  array 
elements.  Consequently,  most  existing  radio  DF  systems  are 
based  on  phase  only  measurements.  As  the  planar  waveform 
impinges  upon  the  receive  elements,  the  difference  in  phase  is 
sampled  (by  several  comparators)  and  indicates  the  direction 
of  arrival.  Ambiguities  in  the  direction  of  arrival  occur  if  the 
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separation  between  interferometer  elements  is  greater  than  half 
the  electrical  wavelength. 

To  resolve  the  ambiguities,  the  addition  of  a  second  interfer¬ 
ometer,  can  be  used  [  1  ] .  Multiple  baseline  systems,  however,  can 
become  quite  long  when  high  resolution  is  required.  Many  ap¬ 
proaches  to  resolving  the  ambiguities  within  a  phase  sampling 
system  have  been  recently  investigated.  The  extended  phase  in¬ 
terferometry  technique  to  resolve  the  ambiguities  incorporates 
both  calibrated  phase  and  amplitude  response  data  from  the  an¬ 
tenna  arrays  [2].  This  is  done  by  appropriately  weighting  the 
square  of  the  baseline  phase  differences  with  the  antenna  gains. 
The  incorporation  of  amplitude  data  provides  significant  perfor¬ 
mance  improvement  over  phase  only  interferometry,  however,  it 
requires  a  modest  increase  in  computational  complexity.  Pseudo 
Doppler  techniques  have  also  been  investigated  in  the  2-2000 
MHz  range  [3].  In  this  approach,  the  system  employs  a  circular 
array  of  four  omni-directional  antennas.  A  commutation  switch 
samples  the  antennas  such  that  the  resulting  phase  modulation 
on  the  received  signal  encodes  the  direction  of  arrival.  As  the 
elements  are  commutated  an  FM  discriminator  responds  with 
impulses  with  size  approximately  proportional  to  the  phase  step 
between  each  element.  The  maximum  spacing  between  adja¬ 
cent  elements  is  less  than  half  a  wavelength.  The  phase  of  the 
output  is  used  to  quantize  an  estimated  angle  of  arrival.  The 
single  channel  interferometer  minimizes  cost  and  power  con¬ 
sumption  by  avoiding  a  multitude  of  expensive  receivers  and  by 
employing  the  cost  effective  pseudo  Doppler  antenna.  The  sym¬ 
metric  pair  antenna  array  provides  accurate  DF  by  means  of  low 
resolution  phase  measurements,  achieved  by  virtue  of  the  sym¬ 
metry  of  the  array  geometry  [4].  This  type  of  array  consists  of 
pairs  of  antenna  elements  that  are  arranged  around  a  common 
phase  center.  Due  to  symmetry,  any  arriving  wavefront  gener¬ 
ates  a  phase  lead  at  one  member  and  a  phase  lag  of  equal  mag¬ 
nitude  at  the  other  member  of  the  symmetric  pair.  Since  the  po¬ 
larity  of  the  output  signal  is  dependent  on  the  angle  of  incidence, 
it  can  be  used  to  determine  the  signal’s  angle  of  arrival.  The 
sign  bit  pattern  of  the  pairs  of  a  given  symmetric  pair  array  will 
change  for  different  angles  of  incidence  and  assuming  an  unam¬ 
biguous  arrangement  (annular  zones  on  the  surface  of  a  sphere) 
must  be  distinguishable. 

This  paper  identifies  a  new  phase  sampling  interferometer  ap¬ 
proach  that  can  provide  a  high  resolution,  small  baseline  array 
with  fewer  number  of  phase  sampling  comparators.  The  ap¬ 
proach  is  based  on  preprocessing  the  received  signal  using  sym¬ 
metrical  number  systems  (SNS).  Antennas  based  on  both  an  op¬ 
timum  symmetrical  number  system  (OSNS)  and  a  robust  sym¬ 
metrical  number  system  (RSNS)  are  investigated.  The  SNS  pre- 
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processing  is  used  to  decompose  the  spatial  filtering  operation 
into  a  number  of  parallel  suboperations  (moduli)  that  are  of 
smaller  computational  complexity.  Each  suboperation  is  a  sepa¬ 
rately  configured  interferometer  that  symmetrically  folds  the  de¬ 
tected  phase  difference  with  folding  period  equal  to  the  twice  the 
modulus  2mi  (for  the  OSNS)  and  2Nmi  (for  the  RSNS)  where 
N  is  the  number  of  interferometers  that  are  used  within  the 
linear  array.  A  small  comparator  ladder  midlevel  quantizes  each 
folded  phase  response.  Consequently,  each  suboperation  only 
requires  a  precision  in  accordance  with  that  modulus.  A  much 
higher  DF  spatial  resolution  is  achieved  after  the  N,  different 
moduli  are  used  and  the  results  of  these  low  precision  subop¬ 
erations  are  recombined.  By  incorporating  the  OSNS  or  RSNS 
preprocessing  concept,  the  field  of  view  of  a  specific  configura¬ 
tion  of  interferometers  and  phase  sampling  comparator  arrange¬ 
ments  can  be  analyzed  exactly.  The  OSNS  gives  the  maximum 
dynamic  range  or  number  of  spatial  resolution  bins  while  the 
RSNS  reduces  considerably  the  number  of  possible  encoding 
errors.  Experimental  results  for  both  a  5-bit  OSNS  and  a  6-bit 
RSNS  array  are  compared.  The  errors  in  the  encoding  of  the  di¬ 
rection  of  arrival  are  quantified  for  both  architectures. 

The  outline  of  this  paper  is  as  follows.  In  Section  II,  two- 
element  interferometry  is  briefly  reviewed.  In  Section  III,  the 
OSNS  is  presented  and  the  practical  implications  of  building 
the  OSNS  DF  antenna  are  addressed.  Experimental  results  are 
shown  to  demonstrate  the  transfer  function  of  the  array  and 
quantify  the  possible  encoding  errors.  Section  IV  presents  the 
RSNS  and  the  design  equations  for  the  RSNS  array.  Experi¬ 
mental  results  are  also  shown  for  comparison.  It  is  shown  that 
although  the  RSNS  has  a  somewhat  smaller  dynamic  range  than 
the  OSNS,  the  absence  of  encoding  errors  makes  it  particularly 
attractive  for  modern  DF  systems. 

II.  Phase  Sampled  Intereerometery  and  Folding 
Waveeorms 

A  two  element  linear  interferometer  is  shown  in  Fig.  1 .  The 
two  antenna  elements  are  spaced  a  distance  d  apart  and  the  in¬ 
cident  plane  wave  arrives  with  bearing  angle  9b- this  phase 
monopulse  configuration  the  angle  is  measured  from  the  per¬ 
pendicular  to  the  baseline  axis  and  can  take  on  values  7r/2  > 
9b  >  — 7r/2.  The  phase  difference  between  the  two  elements  is 

27r 

At/j  —  t/ji  -  tlj2  —  ^  dsm{9B)  (1) 

and  is  a  function  of  the  incidence  angle  of  the  wave. 

The  signals  received  by  the  antenna  elements  are  mixed 
(multiplied  together  and  lowpass  filtered)  resulting  in  an  output 
signal  whose  frequency  is  the  difference  of  the  two  input  signal 
frequencies.  Because  the  signals  have  the  same  frequency,  the 
mixer  output  voltage  is  a  value  proportional  to  the  difference 
between  signal  phases  i/’i  and  i/>2-  This  difference  is  not  purely 
At/j  because  of  the  time  delays  due  to  the  different  transmission 
line  lengths  from  each  antenna  element.  However,  these  line 
lengths  are  known  and  can  be  compensated  for  in  the  angle 
estimate.  Let  the  signals  from  the  two  antenna  elements  be 


(2) 


Fig.  1.  Linear  two  element  interferometer  array  geometry. 


and 


V2{t)  -  V  cos[27r/t  -h  (f>2{t)] 


(3) 


where  V  is  the  maximum  value  of  the  voltage  at  the  antenna  ele¬ 
ments.  Let  i/>o  be  the  phase  difference  between  the  transmission 
lines  to  the  two  elements.  The  low-pass  mixer  output  voltage  is 


f  2'Kd  . 

=  -  cos  -  SIT 


Vout  =  cos(A(/>)  =  —  cos  (  —  sm{9B)  +  V-’o 


(4) 


which  contains  the  plane  wave  angle  of  arrival  (AO A)  informa¬ 
tion.  For  values  of  of  =  A/2  and  t/jQ  —  0,  A(f>  —  7rsin(0B). 
As  the  AOA  9b  varies  from  —ty/2  to  ty/2,  the  phase  differ¬ 
ence  A(/varies  from  — tt  to  tt  as  shown  in  Fig.  2(a).  The  output 
voltage  from  the  phase  detector  is  also  a  function  of  the  phase 
difference  and  is  a  symmetrical  folding  periodic  waveform.  To¬ 
gether  these  relationships  give  the  phase  detector  output  voltage 
as  a  symmetrical  function  of  the  AOA  as  shown  in  Fig.  2(b)  for 
d  =  A/2. 

Ambiguities  are  generated  for  baselines  where  d  >  A/2.  That 
is,  the  phase  detector  output  voltage  is  highly  ambiguous  with 
a  single  phase  corresponding  to  many  angles  of  arrival.  The 
number  of  folding  periods  n  that  occur  within  an  AOA  of  tt  ra¬ 
dians  is 


2d 

T' 


(5) 


For  example,  with  d  —  7.5A,  n  —  15  folds  are  available  as 
shown  in  Fig.  3.  The  folding  period  is  not  constant  but  grows 
larger  in  proportion  to  the  angle  off  broadside  because  of  the 
siii(0B)dependence  in  (4). 

The  ambiguities  within  the  symmetrical  folding  waveforms 
represent  the  phase  difference  between  the  elements  and  can  be 
resolved  by  using  additional  interferometers  in  the  linear  array. 
Typically,  each  interferometer  in  the  linear  array  symmetrically 
folds  the  phase  response  with  the  folding  period  between  in¬ 
terferometers  being  a  successive  factor  of  two  or  ((4  =  2ds  = 
4d2  =  8di  ■  High  speed  binary  comparators  are  used  to  produce 


vi{t)  =  V  cos[27r/f  + 
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Fig.  2.  Two-element  interferometry  with  d  =  A/2 .  (a)  Mixer  output  voltage  versus  phase  difference  between  two  elements,  (b)  Mixer  output  voltage  versus  angle 
of  incidence. 


Fig.  3.  Mixer  output  voltage  versus  AOA  for  d  =  7.5  A. 


a  digital  output.  The  folding  waveforms  are  shifted  appropri¬ 
ately  using  a  phase  shifter  in  each  channel  to  achieve  a  Gray 
code  result.  The  folded  output  from  each  phase  detector  is  then 
quantized  with  a  single  comparator  with  a  normalized  threshold 
level  T  —  0.5.  Together,  the  comparator  outputs  directly  en¬ 
code  the  signal’s  AOA  in  the  Gray  code  format.  This  approach 
makes  use  of  the  periodic  dependence  of  the  interferometer’s 
phase  response  on  the  applied  plane  wave’s  AOA  and  the  dis¬ 
tance  between  the  elements  of  each  interferometer.  One  of  the 


major  limitations  associated  with  this  approach  is  the  achiev¬ 
able  resolution.  For  the  folding  periods  to  be  a  successive  factor 
of  2,  the  distance  between  the  elements  must  also  be  doubled. 
That  is,  an  8-bit  DF  antenna  using  the  previous  scheme  would 
require  element  spacings  A/4,  A/2,  A,  2A,  4A, . . .,  32A  with  a 
total  baseline  length  of  32A.  This  distance  doubling  of  the  ele¬ 
ment  spacings  requires  complex  analog  hardware  and  adversely 
affects  the  physical  implementation  of  the  DF  architecture  and 
ultimately  constrains  the  achievable  resolution. 
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Fig.  4.  Folding  waveforms  and  integer  values  within  moduli  mi  =  5  and  m2  =  6.  (a)  The  OSNS.  (b)  The  RSNS. 


III.  The  OSNS  DF  Antenna 

Incorporation  of  the  OSNS  encoding  overcomes  the  major 
limitations  of  the  binary  encoding  approach  and  provides  an  ef¬ 
ficient  method  to  enhance  the  resolution  of  the  array  while  min¬ 
imizing  the  total  baseline.  The  OSNS  provides  an  efficient  solu¬ 
tion  by  using  additional  baselines  to  correct  ambiguities  in  the 
first  baseline  using  multilevel  quantization  and  prime  sequences 
to  obtain  uniqueness  over  the  dynamic  range. 

A.  Optimum  Symmetrical  Number  System 

The  OSNS  is  composed  of  a  number  of  pairwise  relatively 
prime  (PRP)  moduli  nii.  The  integers  within  each  OSNS  mod¬ 
ulus  are  representative  of  a  symmetrically  folded  waveform  with 
the  period  of  the  waveform  equal  to  twice  the  modulus,  i.e.,  2mi 
[5].  For  m  given,  the  row  vector  below  gives  the  integer  values 
within  twice  the  modulus 

=  [0,  1,  . . . ,  m  -  1,  m  -  1,  . . . ,  1,  0].  (6) 

Fig.  4(a)  shows  part  of  the  OSNS  folding  waveforms  and  the 
integer  values  within  the  modulus  for  nii  —  5  and  m2  —  6. 
The  horizontal  axis  represents  the  normalized  input.  The  vertical 
lines  represent  folding  waveform  reference  levels.  The  numbers 
at  the  top  of  the  figure  represent  the  number  of  reference  levels 


that  are  crossed  by  the  folding  waveform  for  a  given  input  value. 
The  period  of  one  complete  fold  is  equal  to  2mi  =  10  and 
2m2  =  12. 

Due  to  the  presence  of  ambiguities,  the  integers  within  the 
vector  derived  by  (6)  do  not  form  a  complete  system  of  length 
2m  by  themselves.  The  ambiguities  that  arise  within  the  mod¬ 
ulus  are  resolved  by  considering  the  paired  values  from  all  chan¬ 
nels  together.  If  N  pairwise  relatively  prime  moduli  are  used, 
the  dynamic  range  (the  number  of  unambiguous  vectors)  of  this 
scheme  is  [6] 

N 

M  (7) 

i=l 

This  dynamic  range  is  also  the  position  of  the  first  repetitive 
moduli  vector.  As  shown  in  Fig.  4(a),  the  dynamic  range  for  the 
mi  =  5  and  m2  =  6  case  is  M  —  30.  That  is,  no  set  appears 
twice  for  a  normalized  input  in  the  range  [0,  29]. 

B.  OSNS  Antenna  Architecture 

Fig.  5  shows  the  schematic  diagram  of  an  =  2  channel 
OSNS  DF  antenna.  The  array  consists  of  three  elements  with 
all  channels  sharing  a  common  element.  With  the  proper  dis¬ 
tances  between  elements,  each  interferometer  folds  the  phase 
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Fig.  5.  Schematic  diagram  of  the  OSNS  array  with  mi  =  6  and  m2  =  11. 


response  at  2mi.  The  phase  response  from  each  channel  is  quan¬ 
tized  using  nii  —  1  comparators.  The  integers  shown  on  top  of 
both  folding  waveforms  in  Fig.  4(a)  represent  the  number  of 
comparators  ON  due  to  the  phase  voltage  exceeding  the  com¬ 
parator  matching  threshold  voltage. 

The  distance  between  each  pair  of  elements  is  a  function  of 
the  modulus,  m^.  The  phase  response  must  have  the  correct 
number  of  folds  to  cover  the  dynamic  range.  For  a  given  mod¬ 
ulus  nii,  the  number  of  folds  within  the  dynamic  range  is 

rii  =  M/2mi.  (8) 


Since  each  fold  of  the  array  corresponds  to  a  change  in  phase 
of  TT  for  d  —  A/2,  the  required  spacing  between  the  reference 
element  and  the  respective  element  of  channel  i  is  given  by 


MX 

Ami 


(9) 


The  bin  widths  are  uniform  in  size  in  sin(0B)  space  with  a  width 
2/M.  Therefore,  the  spatial  resolution  is  given  by 


Tk  =  arcsm 


2k -M  +  2 


2k -M 


where  A:  £  {0,  1,  ...  M  —  1}  and  A;  =  0  is  the  bin  that  starts  at 
9b  —  — 7r/2.  The  variable  ^  is  a  scale  factor  that  can  be  used  to 


compress  the  encoded  field  of  view  symmetrically  about  broad¬ 
side.  The  phase  waveforms  are  not  ideal  over  180°  because  of 
the  element  pattern  degradation  at  wide  angles.  Since  the  wave¬ 
forms  at  wide  angles  are  not  usable,  the  SNS  can  be  remapped  to 
the  usable  portion  of  the  waveform  to  increase  resolution.  The 
scale  factor  is  defined  to  be 


d!_  _  sin(6') 
d  sin(6'') 


(11) 


where  9  corresponds  to  phase  angle  for  the  unsealed  configura¬ 
tion  and  9'  is  the  desired  phase  angle  for  the  compressed  config¬ 
uration.  Note  that  the  spatial  resolution  is  not  uniform  over  the 
range  of  — 7r/2  to  7r/2.  The  resolution  is  finer  than  r  at  broad¬ 
side  {9  —  0°)  and  increases  with  angle  off  broadside,  because 
of  the  sin(0B)  dependence  in  (5). 


C.  Experimental  Results 

A  linear  array  based  on  the  moduli  mi  =  6  and  m2  =  11 
was  designed,  fabricated,  and  tested  at  a  frequency  8.5  GHz. 
The  radiating  elements  are  open-ended  waveguides  as  shown  in 
Fig.  5.  For  moduli  mi  =  6  and  m2  =  11,  di  —  9.70  cm  and 
d2  —  5.28  cm.  With  M  —  66,  the  number  of  folding  periods 
m  —  5.5  and  712  =  3.  To  provide  an  adequate  signal-to-noise 
ratio,  a  low-noise  amplifier  is  included  at  the  output  of  each 
interferometer  element.  Since  the  common  element  splits  the 
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Channels  are  phased  up  at  50  degrees  where  OSNS 
code  starts  (right-to-left) 


Fig.  6.  Measured  mixer  output  voltages  for  the  OSNS  array. 


signal  into  N  paths,  an  attenuator  is  placed  in  the  other  branches 
to  balance  the  amplitudes.  The  amplifiers  operate  in  saturation 
so  that  the  mixer  input  signal  levels  are  independent  of  the  angle 
of  incidence.  A  fixed  phase  shifter  is  also  included  in  one  branch 
of  each  interferometer  so  that  the  symmetrically  folded  phase  re¬ 
sponse  waveforms  from  each  mixer  may  be  aligned.  This  align¬ 
ment  insures  that  the  comparators  in  the  digital  processor  prop¬ 
erly  sample  the  phase  waveform  and  encode  it  in  the  OSNS. 

The  anechoic  chamber  facility  at  the  Naval  Postgraduate 
School  was  used  for  the  antenna  pattern  measurements  for  each 
pair  of  elements.  The  measured  mixer  outputs  are  shown  in 
Fig.  6.  Since  the  array  is  not  useful  in  the  endfire  regions,  the 
phase  response  from  each  channel  was  minimized  at  an  input 
AOA  of  —50°.  This  is  the  point  where  the  digital  code  starts. 
Introducing  the  appropriate  phase  shifts  at  a  point  between  the 
elements  and  the  mixer  controls  the  locations  of  the  voltage 
minima.  A  simulation  of  the  transfer  function  for  the  OSNS 
array  is  shown  in  Fig.  7(a)  for  ^  =  1.  The  measured  transfer 
function  (output  of  digital  hardware)  is  shown  in  Fig.  7(b). 

The  spikes  in  the  transfer  function  are  AOA  encoding  errors. 
The  phase  waveforms  must  cross  the  comparator  thresholds 
(code  transitions)  at  the  same  time  for  all  channels.  The  OSNS 
can  present  a  problem  at  each  of  these  specific  AOAs.  When 
some  comparators,  at  positions  to  change,  do  change,  while 
others  do  not,  the  recovered  amplitude  will  have  a  large  error. 
Encoding  errors  can  also  occur  due  to  imperfections  in  the 
transmission  line  phases.  These  errors  however,  can  be  isolated 
with  additional  comparators  and  corrected  using  interpolation 
circuits  as  described  in  [6].  An  alternate  approach,  however,  is 
to  develop  a  symmetrical  number  system  that  eliminates  the 
possibility  of  this  type  of  encoding  error.  That  is,  a  symmetrical 


number  system  where  only  one  comparator  threshold  is  crossed 
at  any  one  particular  code  transition. 

IV.  The  RSNS  Antenna 

The  RSNS  preprocessing  architecture  is  a  modular  scheme 
in  which  the  integer  values  within  each  modulus  (comparator 
states),  when  considered  together,  change  one  at  a  time  at  the 
next  position  (Gray  code  properties)  [7].  Although  the  observed 
dynamic  range  of  the  RSNS  is  less  than  that  of  the  OSNS  given 
in  (7),  the  RSNS  Gray  code  properties  make  it  particularly  at¬ 
tractive  for  error  control.  With  the  RSNS  preprocessing,  the  en¬ 
coding  errors  due  to  comparator  thresholds  not  being  crossed 
simultaneously  are  eliminated.  As  a  result,  the  interpolation  cir¬ 
cuits  can  be  removed  and  only  a  small  number  of  comparators 
are  required. 

A.  The  Robust  Symmetrical  Number  System 

In  the  RSNS,  N  different  periodic  symmetrical  wave¬ 
forms  are  used  with  pairwise  relatively  prime  integers 
rtii,  m2,  ■  ■  ■ ,  mj\T.  The  RSNS  is  based  on  the  following 
sequence: 

x„^l  -  [0,  1,  2,  . . . ,  m  -  1,  m,  m  -  1,  . . . ,  2,  1]  (12) 

where  x„ii  is  ^  row  vector  and  m  is  a  positive  integer  (m  >  0). 
In  an  iV-channel  RSNS,  where  N  >  2,  the  basic  sequence  for 
the  (ith  channel  (modulus  m)  is 

Xml  =  [0,  0,  . . . ,  0,  0,  1,  1,  . . . ,  1,  1,  . . . ,  m 

m,  . . . ,  m,  m,  . . . ,  1,  1,  . . . ,  1,  1].  (13) 
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Fig.  7.  Transfer  function  for  the  OSNS  array,  (a)  Simulation,  (b)  Measured  results. 


In  this  sequence  each  value  in  the  x„ii  row  vector  is  put  N 
times  in  succession.  This  sequence  is  repeated  in  both  direc¬ 
tions,  forming  a  periodic  sequence  with  the  period 

■Prsns  =  2mN.  (14) 

Considering  a  single  channel,  the  discrete  states  of  the  robust 
symmetrical  number  system  can  be  expressed  as 

^  ■ 

^  ,  Si  <n<  Nnii  +  Si  +  1 

2Nmi  +  N-  n  +  Si-l 

iV  _ 

Nnii  +  Si  +  2  <  n  <  2Nmi  +  Si  —  1 

(15) 


where 

g  nth  term  of  channel  ?; 

mi  channel  modulus; 

Si  corresponding  sequence  shift  Si  =  0,  1,  2,  N  —1 
(mod  N)  and  TV  >  2  is  the  number  of  channels  in  the 
system. 

The  values  {si,  S2,  . . . ,  sj\t}  must  form  a  complete  residue 
system  modulo  N.  Because  of  the  relative  property  of  the  shifts, 
one  of  the  shift  values  will  be  set  equal  to  0.  The  index  n  cor¬ 
responds  to  the  input  value.  The  discrete  states  of  the  RSNS  are 
indexed  using  the  Si  =  0  row  vector  with  the  index  starting  from 
the  first  zero.  Note  that  the  largest  integer  within  each  periodic 
sequence  is  the  modulus  m. 

An  iV-channel  RSNS  is  formed  of  vectors  by  picking  N 
moduli  mi  and  N  shift  values  Si,  1  <  i  <  N.  Since  the 
fundamental  period  for  channel  i  is  2Nmi,  it  follows  that  the 
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period  for  the  RSNS  vectors  must  be  a  multiple  of  2Nmi. 
Therefore,  the  fundamental  period  for  the  RSNS  is 

^Aisns  =  [2mi7V,  2m2N,  ■■■,  2m ^N]  (16) 

where  [ai,  (12 ^  aA^]  is  the  least  common  multiple  of 
ai,  (12 ^  ,  (IN-  From  number  theory  [7] 

F’i^RSNS  =  2iV[mi,  m2,  . . . ,  mAr],  (17) 


different.  When  the  channel  moduli  are  of  the  form  mi  =  2^  —  1, 
m2  =  2^,  m3  =  2^  +  1,  the  closed-form  expression  for  Mrsns 
is  conjectured  to  be 

4Frsns  =  ftRi  +  ^mi  -h  7  (18) 

where  mi  >  3.  For  two  channels,  with  the  moduli  separated  by 
3  or  more 

Mrsns  =  4mi -I- 2m2  -  2.  (19) 


To  illustrate  the  RSNS,  Fig.  4(b)  shows  the  folding  waveform 
and  integer  values  within  each  modulus  for  mi  =  5  (shift 
Si  =  0)  m2  =  6  (shift  S2  =  1).  Note  the  Gray-code  properties 
from  one  code  transition  to  the  next.  The  thresholds  shown  on 
the  vertical  axis  represent  the  integer  values  within  each  RSNS 
modulus.  The  integer  values  occur  N  =  2  times  in  succession. 

The  system  dynamic  range  (Mrsns)  of  the  RSNS  is  the  max¬ 
imum  number  of  distinct  vectors  without  a  redundancy.  The 
performance  of  the  2-channel  example  is  shown  in  Fig.  4(b). 
The  selection  of  the  shift  and  the  permutations  between  the  2 
moduli  has  no  effect  on  Mrsns  -  However,  the  point  indices  cor¬ 
responding  to  Mrsns  (beginning  point  and  ending  point)  are 


Checking  against  computer  results  when  mi  =  8  and  m2  =  17, 
M  =  64  which  fits  (19). 


B.  RSNS  Antenna  Design 

To  demonstrate  the  efficiency  of  the  RSNS  preprocessing,  the 
design  of  a  6-bit  antenna  using  mi  =  8  and  m2  =  17  is  con¬ 
sidered.  A  schematic  diagram  of  the  RSNS  antenna  is  shown  in 
Fig.  8  for  5i  =  0,  52  =  1  resolution  bins.  The  number  of  folds 
within  each  modulus  is 


M  _  M 
-Prsns  2miN 


(20) 
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Fig.  9.  RSNS  logic  block. 


or  ni  —  2,  and  n2  —  0.94.  The  required  distance  between 
antenna  elements  is 


angles  a  reporting  error  exists.  The  AOA  resolution  or  bin  width 
in  9-space  for  the  kth  bin  is 


MX 

4miN 


(21) 


Tk  =  arcsm 


2k  -  M  +  2 


—  arcsm 


2k -M 

.  CM  , 


(23) 


With  a  scale  factor  of  ^  =  2/a/3,  di  —  4.33  cm,  and  d2  — 
2.03  cm.  For  a  normalized  folding  waveform  amplitude  (— 1 
to  1),  the  threshold  for  the  jth  comparator  for  the  modulus  mi 
channel  is 

{  mi  —  j  -\-  h  \ 

Vj^rrii  =  cos  (  - ^  TT  )  .  (22) 

\  rrii  J 

The  input  signal  is  applied  in  parallel  to  both  interferometers. 
The  output  from  each  phase  detecting  mixer  is  amplitude 
analyzed  using  mi  comparators.  For  this  design,  there  are 
a  total  of  25  comparators  with  a  maximum  of  17  loaded  in 
parallel.  An  RSNS-to-binary  logic  block  translates  the  RSNS 
residues  (number  of  comparators  ON  in  each  thermometer 
code)  to  a  more  convenient  representation.  The  logic  block  and 
corresponding  index  n  are  shown  in  Fig.  9. 

A  consequence  of  quantizing  the  angle  of  arrival  into  a  bin  is 
that  the  RSNS  system  reports  any  signal  that  falls  within  a  bin 
as  if  it  arrived  at  the  bin  center,  9k-  That  is  to  say,  for  one  angle 
within  each  bin,  the  estimation  is  exact,  but  for  the  remaining 


Another  advantage  of  the  RSNS  is  that  small  phase  errors  can 
be  tolerated  without  a  serious  error  to  a  reported  AOA.  These 
phase  errors  cause  the  folded  waveform  from  the  mixer  to  shift, 
thus,  changing  the  timing  of  the  comparator  state  changes.  Be¬ 
cause  the  waveforms  from  the  mixers  are  quantized,  the  bin 
width  or  Vk  provides  some  freeplay  in  the  antenna  circuitry  and 
limits  AOA  reporting  errors  to  the  adjacent  quantization  levels 
(bins).  However,  there  is  a  limit  to  this  feature.  If  the  phase  error 
becomes  too  great,  the  comparators  will  not  change  state  in  the 
correct  sequence,  which  may  or  may  not  correspond  to  a  code 
that  is  mapped  in  the  dynamic  range.  These  errors  are  similar 
to  the  errors  that  the  OSNS  experiences  when  the  comparators 
of  both  channels  do  not  change  state  at  the  same  time.  For  the 
system  designed,  a  theoretical  phase  error  of  6°  can  be  toler¬ 
ated  in  the  modulus  17  channel  and  a  phase  error  of  12°  in  the 
modulus  8  channel.  Figs.  10  and  11(a)  and  (b)  show  the  simu¬ 
lated  transfer  function.  Fig.  10  is  the  transfer  function  with  no 
phase  error  in  both  channels.  Fig.  1 1(a)  and  (b)  are  the  modulus 
8  channel  with  zero  degrees  of  phase  error  and  3°  and  6°  phase 
error  in  the  modulus  17  channel.  Note  how  the  transfer  function 
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Fig.  10.  Simulated  transfer  function  with  no  phase  errors  in  either  channel. 


Fig.  1 1 .  Modulus  17  channel,  (a)  3“^ .  (b)  6“^  phase  error  introduced. 

becomes  jagged  and  erratic  as  the  phase  error  is  increased  until 
it  becomes  discontinuous  at  6°. 


C.  Experimental  Results 

A  RSNS  array  based  on  the  moduli  mi  =  8  and  m2  —  11  was 
designed,  fabricated,  and  tested  at  a  frequency  of  8.0  GHz.  The 
schematic  diagram  of  the  two-channel  array  is  shown  in  Fig.  8. 
The  radiating  elements  are  printed  circuit  dipoles.  For  moduli 
mi  =  8  and  m2  =  17,  the  unsealed  distances,  ^  =  1,  are 
di  —  3.76  cm  and  d2  —  1.78  cm.  With  M  —  64,  the  number 
of  folding  periods  rii  —  2  and  n2  —  0.94.  To  provide  an  ad¬ 
equate  signal-to-noise  ratio,  a  low-noise  amplifier  is  included 
at  the  output  of  each  interferometer  element.  Also  included  in 
the  RSNS  array  is  a  lowpass  filter  in  each  channel  to  eliminate 
the  amplifier  harmonics  that  are  present.  Since  the  common  ele¬ 
ment  splits  the  signal  into  N  paths,  an  attenuator  is  placed  in  the 
other  branches  to  balance  the  amplitudes.  A  fixed-phase  shifter 
is  also  included  in  one  branch  of  each  interferometer  so  that 
the  symmetrically  folded  phase  response  waveforms  from  each 
mixer  may  be  aligned.  This  alignment  insures  that  the  compara¬ 
tors  in  the  digital  processor  properly  sample  the  phase  wave¬ 
form  and  encode  it  in  the  RSNS.  For  the  mi  =  8  channel,  8 
comparators  are  required.  For  the  m2  =  17  channel  17  com¬ 
parators  are  required.  The  electrically  erasable  read  only  mem- 
bory  (EEPROM)  that  was  used  accommodates  15  inputs.  The 
modulus  17  channel  processing  is  conducted  in  two  stages.  In 
the  first  stage,  the  4-bit  thermometer  to  binary  encoder  maps 
the  15  least  significant  codes  out  of  the  17  comparators  (ones 
with  the  smallest  threshold).  The  remaining  two  comparators 
are  mapped  on  the  second  EEPROM  along  with  the  4-bit  binary 
encoder  output  and  all  modulus  8  comparators.  The  digital  pro¬ 
cessor  puts  out  a  bin  number  k  corresponding  to  the  estimated 
AOA  9.  The  estimated  AOA  for  bin  k  is  given  as 


9k  -  sin  ^ 


k2k  +  l 

V  m 


(24) 


Broadside  incidence  for  the  RSNS  antenna  corresponds  to  the 
transition  between  Bins  3 1  and  32.  Eor  the  modulus  17  channel, 
at  slightly  negative  angles  (bin  31),  the  ninth  comparator  is  ON 
and  for  slightly  positive  angles  (bin  32),  the  ninth  comparator 
is  OPE.  That  is,  the  phase  value  at  9b  —  0  corresponds  exactly 
to  the  threshold  of  the  ninth  comparator.  Prom  (22)  the  phase 
corresponding  to  the  threshold  value  F17, 9  is  tt/2.  This  is  the 
phase  value  required  for  the  modulus  17  channel  at  9b  =  0 
and  is  obtained  using  the  phase  adjuster  shown  in  Pig.  8.  In  the 
modulus  8  channel,  9b  —  0  occurs  midway  between  comparator 
7  and  6  thresholds.  Since  Vs, 7  —  cos“^(57r/16)  and  = 
cos“^(37r/16),  the  required  phase  value  is  47r/16  =  7r/4and  is 
also  obtained  using  the  modulus  8  phase  adjuster. 

The  normalized  mixer  folding  waveform  outputs  are  shown 
with  the  simulated  waveforms  in  Pig.  12(a)  for  the  mi  =  8 
channel  and  Pig.  12(b)  for  m2  =  17  channel  for  ^  =  1.0  and 
the  basic  features  correspond  to  predicted  curves.  The  measured 
folding  outputs  from  the  shift  and  bias  amplifiers  in  both  chan¬ 
nels  are  shown  in  Pig.  13.  The  digital  circuit  maps  the  phase 
difference  at  the  mixer  output  into  an  estimated  angle  of  arrival. 
The  predicted  output  of  the  unsealed  prototype  array  is  repro¬ 
duced  as  Pig.  14(a).  The  measured  mixer  output  contains  phase 
errors  that  result  in  differences  between  the  simulated  and  pre¬ 
dicted  data.  These  phase  errors  cause  the  thermometer  code  to 
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Fig.  12.  Normalized  mixer  output  and  simulation  waveforms  for  RSNS  array, 
(a)  m.i  =  8.  (b)  m2  =  17. 

incorrectly  map  the  incident  wave  to  the  proper  angle  of  arrival. 
The  measured  transfer  function  of  the  unsealed  prototype  an¬ 
tenna  is  shown  in  Fig.  14(b). 

The  large  discontinuities  in  the  transfer  function  are  due  to 
phase  and  amplitude  errors  that  arise  from  several  sources.  Mu¬ 
tual  coupling  is  a  common  cause  of  phase  front  distortion  in 
small  arrays  of  closely  spaced  elements.  In  order  to  increase  the 
distance  between  elements  and  thus  reduce  the  mutual  coupling, 
a  scale  factor  ^  can  be  introduced.  The  scale  factor  is  the  ratio 
of  the  element  spacings  of  the  scaled  and  unsealed  arrays.  Re¬ 
ducing  the  scale  factor  narrows  the  mapable  field  of  view  of  the 
antenna  and  increases  the  resolution  within  the  mapable  field  of 
view.  For  a  scale  factor  of  ^  =  2/ a/3,  the  maximum  mapable 
aperture  (Omm)  is  60°.  The  predicted  transfer  function  for  the 
scaled  array  is  shown  in  Fig.  15(a).  The  corresponding  measured 
transfer  function  is  reproduced  as  Fig.  15(b).  The  mixer  outputs 
contain  significant  phase  errors,  and  the  increased  resolution  of 
the  antenna  limits  the  ability  to  process  the  errors.  The  steeper 


Fig.  13.  Measured  folding  waveform  outputs  from  the  shift  and  bias 
amplifiers. 


(a) 


(b) 

Fig.  14.  RSNS  transfer  functions,  (a)  Predicted,  (b)  Measured. 
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Fig.  15.  RSNS  transfer  functions,  scale  factor  {  =  2/^/3.  (a)  Predicted, 
(b)  Measured. 

slope  of  the  scaled  transfer  function  is  clearly  visible  when  com¬ 
paring  Figs.  14(a)  and  15(a). 

The  transfer  function  for  the  scaled  array  of  dipoles  does  not 
show  any  significant  improvement  due  to  the  reduction  of  mu¬ 
tual  coupling.  To  verify  this,  a  second  array  comprised  of  open 
ended  waveguide  elements  was  constructed  and  tested.  These 
waveguide  elements  have  significantly  less  H-plane  mutual  cou¬ 
pling  than  the  dipoles,  but  the  transfer  function  showed  no  im¬ 
provement  over  with  printed  circuit  (dipole)  elements.  This  indi¬ 
cates  that  mutual  coupling  is  not  the  primary  cause  of  the  phase 
errors.  The  microwave  circuit  has  three  likely  sources  of  errors: 
1)  The  low  noise  amplifiers;  2)  coaxial  cables  and  phase  ad¬ 
justers;  and  3)  the  rigid  connectors.  The  prototype  uses  two  cas¬ 
caded  LNAs  operating  in  saturation  in  each  channel.  The  second 
amplifier  amplifies  harmonics  generated  by  the  first  amplifier. 
The  second  amplifier  also  introduces  its  own  harmonics.  Filters 
were  used  to  suppress  harmonics. 


The  coaxial  cables  and  phase  adjusters  are  another  source  of 
error  within  microwave  circuits.  The  rigid  coaxial  cables  are 
coarsely  trimmed  at  8.0  GHz  and  the  phase  shifters  are  used  for 
fine  trimming.  The  cable  lengths  are  long  and  subject  to  thermal 
expansion  and  contraction,  particularly  in  the  vicinity  of  the 
amplifiers.  Finally,  all  of  the  devices  in  the  circuit  are  joined  with 
SMA  type  connectors.  The  circuit  was  repeatedly  assembled  and 
disassembled  for  troubleshooting  and  consequently  the  quality 
of  the  connections  was  degraded.  Phase  errors  of  ±2°  can  easily 
be  introduced  at  each  connection.  Thus,  it  appears  that  many 
small  errors  accumulate  in  a  manner  such  that  the  6°  and  12°  error 
thresholds  are  exceeded,  resulting  in  a  degraded  transfer  function. 

V.  Conclusion 

This  research  has  introduced  the  OSNS  and  the  RSNS  and 
shown  they  can  be  used  to  resolve  ambiguities  in  phase  sampling 
interferometry.  A  direction  finding  antenna  architecture  based  on 
the  OSNS  has  many  inherent  advantages  over  conventional  di¬ 
rection  finding  techniques  including:  simple  microwave  beam¬ 
forming  network  that  incorporates  wideband  components,  wide 
instantaneous  field  of  view,  and  high  resolution  angle  of  arrival 
with  a  variety  of  element  arrangements  having  very  short  base¬ 
lines.  The  design  equations  for  comparator  thresholds,  broadside 
phase  difference,  and  a  scale  factor  are  derived.  The  simulated 
transfer  functions  demonstrate  that  the  OSNS  and  RSNS  antenna 
will  have  a  wide  instantaneous  field  of  view  and  high  resolution. 

Test  arrays  were  built  and  their  performance  measured.  The 
measured  transfer  functions,  while  exhibiting  the  basic  features 
of  the  ideal  transfer  function,  contained  significant  encoding  er¬ 
rors.  The  errors  are  attributed  to  phase  errors  in  the  microwave 
circuit  and  its  components,  as  well  as  mutual  coupling.  Using 
a  printed  circuit  design  with  surface  mount  components  can  re¬ 
duce  these  errors  significantly. 
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